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Abstract 

In this paper we compute exactly the ground state energy and entropy of the dilute ferromagnetic 
Ising model. The two thermodynamic quantities are also computed when a magnetic field with 

random locations is present. The result is reached in the replica approach frame by a class of 

i i 

replica order parameters introduced by Monasson [1]. The strategy is first illustrated considering 
the SK model, for which we will show the complete equivalence with the standard replica approach. 
Then, we apply to the diluted ferromagnetic Ising model with a random located magnetic field, 
which is mapped into a Potts model. This formalism can be, in principle, applied to all random 
systems, and we believe that it could be of help in many other contexts. 
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1 - INTRODUCTION 



The theoretical modehng of statistical systems in many areas of physics makes use of 
random Hamiltonians. Assuming self-averaging, observable quantities may be evaluated 
by the logarithmic average of the partition function <log(Z)> over the quenched random 
variables. In this way, the technically difficult task of computing Z for a given realization 
of the disorder is avoided. However, the mathematical operation of directly computing 
<log(Z)> is also difficult, and can be done only using replica trick, which implies the 
computation of <Z">. 

Despite its highly successful application to the treatment of some disordered systems (the 
most celebrated is the SK model [2, Q], the replica trick encounters serious difficulties when 
applied to many other models. The main reason is the unbounded proliferation of replica 
order parameters, as for example the multi-overlaps in dilute spin glass ^, 0] . 

In this paper we use a class of replica order parameters firstly introduced by Monasson [l| , 
which, in principle can be used for all models. We preliminarily observe that the computation 

of <Z"> implies the sum over all the realizations of the spin variables cTj = {o\^ af, , a^), 

any of them may take 2" values. Assume that x{a)N is the number of vectors which 

equal the given vector a = (cr^,cr^, jO"""), then, <Z"> can be re-expressed in terms of 

a sum over the possible positive values of the 2" magnetization densities x(cr), with the 
constraint I]o-x(cr) = 1. In practice, this is equivalent to maximize with respect to these 
order parameters. 

The strategy proposed here will be illustrated, considering the SK model, in next section, 
where we will also show the complete equivalence with the standard approach. Nevertheless, 
all models can be described in terms of the 2" order parameters x{a). We will tackle the 
dilute ferromagnetic Ising model in section 3 and in section 4 we will compute exactly the 
ground state energy and entropy. The results coincide with those found in [6] where the 
approach is not based on replicas. In section 5 we extend the scope by considering the same 
model in presence of a magnetic field with random locations, and, finally, in section 6 we 
are able to compute the ground state energy and entropy also in this case. Conclusions and 
outlook are in the final section. 
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2 - SK SPIN GLASS 



In this section we try to illustrate our approach considering the SK model. We do not have 
new results concerning SK, but we just show the complete equivalence with the standard 
replica approach. 

The partition function is 




where the sum J2# goes over the 2^ realizations of the cTj, the sum J2i>j goes over all 
N{N — l)/2 pairs ij and the Jij are independent random variables with mean and variance 
1. Then, by replica approach, neglecting terms which vanish in the thermodynamic limit, 
we have 

<Z«>=Eexpf|^E(E^r^")0 (2) 

where the sum X)# goes over the 2"^ realizations of the af and the sum J^ij goes over all 
iV^ possible values of i and j. 

Assume that Nx{cr) is the number of vectors {al, af, , cr[^) which equal the given vector 

a = ((T^, cr^, , o""), then, according to we can write 



i.j \ a / (T,T 



x{cr)x{T){aT)'^ (3) 



where J2a,T goes over the 2^" possible values of the variables a and r and where ar is the 
scalar product ar = J2a cr"'T°'- Then, observe that the number of realizations corresponding 
to a given value of the 2^ magnetization densities x{a) is 

exp(-A^ log(x((T)) (4) 

where J2a is the sum over the 2" possible values of the variable a. Indeed, in the above 
expression we neglected terms which are in-influent in the thermodynamic limit. 
We can now define as the large limit of ;^log <Z">, then 

= max [ y ^x((T)x(r)((Tr)2 -^a;((T)log(a;(a))] (5) 
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The maximum is taken over the possible values of the 2" order parameters x{a) provided 
that J2a^{'^) = 1- The constraint can be accounted by adding the Lagrangian multiplier 
\(J2a^i(^) ~ 1) to expression then the maximum is given by 

yE^MM'-M^(^))-i = a (6) 

where A has to be chosen in order to have J2a = 1- From this equation we get that the 
maximum is realized for the set of the 2" order parameters x{a) which satisfy 




where the sum J2t goes over the 2" possible values of the variable r and where A is 

^ = Eexp(f E^(^)K)') (8) 

and where X^o- is the sum over the 2" possible values of the variable a. 

The explicit solution could be found, in principle, by a proper choice of the parametriza- 
tion of the x{a), but it can be easily seen that this solution coincides with the standard 
solution of the SK model. In fact, if one defines the symmetric matrix qajj as 



q^p = Y.x{rKT^ (9) 

T 

one gets 



A 

where 



x(a) = iexp\f3^ ^ | (10) 

a>/3 



i = Eexp E?«/5^"^^ (11) 

o- \ a>l3 J 

the sum J2a>i3 goes over all n{n — l)/2 pairs a, [3 and the diagonal therms of the matrix 
disappeared since they cancel out in the expressions ([7]), ([8]). 
Then, after some work, the expression ([5]) rewrites as 

= max [ ^ n - ^ ^ {q^pf + log ( J] exp(/3^ Q^^^"^^) \ ] (12) 
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where the maximum is over the variables qafs- This is the standard solution of [2] provided 
the proper maximum of qap is found [3.]- Then $ = lim„^o ^ — ^ ~ where 5* is the 
entropy and E the energy. 



3 - DILUTE FERROMAGNET 



In this section, we show how our approach works for the dilute ferromagnetic system. 
This model is much less studied than SK, and many informations about its phenomenology 
are still missing. 

The partition function is 



Z = ^ exp {(3Y^ KijOiUj (13) 
# V i>i ) 

where Kij are quenched variables which take the value 1 with probability and otherwise. 
The dilution coefficient 7 may take any positive value and the number of bonds is about 
This model has been recently studied in a recent paper [5], while models with lesser 



dilution 
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number of bonds of order A^*^ with e > 1 ) have been also considered. 



We can rewrite the above expression as 



Z = exp log(cosh(/3))) E 11(1 + tanh(/5)ir,,or,a,) (14) 

where the equality holds in the sense that \og[Z) coincide in the thermodynamic limit for 
([T3D and ([HD because Ei>j ^ij = + o{N) 
Let us define 



^ = E + tanh(/5)K,,a,a,) (15) 



than we can compute <G"> 



< >= E n f 1 - ^ + i 11(1 + tanh(/?)<a;)) 



(16) 



which is thermodynamically equivalent to 



< >= exp(-^) E ^ E n (1 + tanh(/?)ara^) (17) 



# i>3 
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By introducing the x{a) with identical meaning as in previous section, and defining 
as the large N limit of log <G"> we get 



vl/„ = max [ -| + I E ^(^)^(^) + tanh(/3)a°r°) - ^ ^^(a) log(a;(a)) ] (18) 

(T,r a a 

where the maximum is taken with respect the 2" order parameter x{a). To obtain an explicit 
expression one should find their parametric expression. This will be done in next section 
only for the zero temperature case. Formally, the maximum is reached for 



a) = jexp (^7E^Mn(l + tanh(/?)a"r")j (19) 



where 



A = E exp ( 7 E ^i^) 11(1 + tanh(/3)a"r^ 
Then, according to ([H]), ([I5]), (1171) and ([IH]), we get 



(20) 



$ = lim — + ^log(cosh(/3)) = 5-/3E (21) 
"^0 n 2 

where S is the entropy and E the energy. 

4 - DILUTE FERROMAGNET: ZERO TEMPERATURE 

Let us consider the simpler case of vanishing temperature. In this limit tanh(/3) = 1 and 
expression f|T8|) becomes 

= max [ + ^2" E - E ^(^) log(^(^)) ] (22) 

which is a standard 2"-components Potts model. The solution is known and can be found 
assuming that 2" — 1 quantities x{a) take the value and one takes the value ^-'^^^tt^- 
The state with different value can be any of the possible 2", we assume that is the one with 
a" = 1 for all a. We can write: 

^ ' 2" 2" 



6 



which obviously satisfy the constraint J2a^i'^) = 1- Inserting the above expression in (122|) 
we obtain 

^„ = max[J(2"-l)^^2 + 5„] (24) 

with 

Bn = - (^^^^!^) log(l + (2" - 1)^^) - (2" - log(l -e) + n log(2) (25) 

If we expand to the first order in n the above expression, we obtain 

^„ = max [nlog(2)-^2^nlog(2)(l-0)(l-log(l-0))] (26) 
9 2 

Then, if we take into account fl?I]) and we also take into account that for large (3 one has 
log(cosh(/5)) = P — log(2), we have that the energy at temperature equals 

E = -l (27) 
while the entropy S* is — | log(2) + lim„^o Therefore: 

S = log(2)|(e,2 _ ^) ^ i^g^2)(l - e,){l - log(l - e,)) (28) 
where 9c is given by the equation 

exp(-7^^,) = l-9c (29) 

This equation has a single solution 6'c = if 7 < 1 and one more non trivial solution if 7 > 1 
which corresponds to the maximum. Therefore, at temperature, for 7 < 1 the system is in 
a paramagnetic phase while for 7 > 1 is in a disordered ferromagnetic phase. The transition 
corresponds to the percolation transition generated by the ferromagnetic links. The entropy 
S and the order parameter 6c are plotted in Fig. 1 as a function of the dilution coefficient 7. 
At the transition value 7 = 1, the first derivatives of both entropy and 6c are discontinuous. 

5 - DILUTE FERROMAGNET IN A MAGNETIC FIELD 

In this section, we show how our approach works for the dilute ferromagnetic system in 
a magnetic field with random locations. 
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FIG. 1: Entropy (dashed line) and order parameter 6c (full line) as a function of the dilution 
coefficient 7 at temperature and magnetic field. The transition is at 7 = 1 where the first 
derivatives of both entropy and 6c are discontinuous. 

The partition function of this model is 

Z = ^ exp /5 ^ KijaiCj + /3 ^ hiai (30) 

# V i>j i J 

where Kij are the previously defined quenched variables and the hi take the value h with 
probability 5 and otherwise. We can rewrite the above expression as Z = PG where 



(31) 



and where 



p = exp log(cosh(/5))^ exp (5A^log(cosh(/3/i))) (32) 

the equality holds in the sense that \og{Z) coincide in the thermodynamic limit when 
computed from (!30l) and from Z = PG with G and P given by (|3T|) and (!32|) . 
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If one takes into account that 



< n (1 + tanh(/3/ii)af ) >=l[ll-6 + 6l[{l + tanh(/5/i)af ) j 
one has that \E'„, with respect to (fT8|) . contains the extra term 



(33) 



^x(a)log ll - 6 + 6YI{1 + tanh{(3h)a^)] (34) 

a \ a J 



which is associated to the magnetic field. 

Then, if the correct maximum is found we get 



$ = hm — + - log(cosh(/5)) + 6 log(cosh(/3)) = S - pE (35) 



6 - DILUTE FERROMAGNET IN A MAGNETIC FIELD: ZERO TEMPERATURE 

In the vanishing temperature hmit one has tanh(/5) = tanh(/5/i) = 1 and expression 
becomes 



= max [- 1 + ^2" + E ^(^) log(l - 5 + 5 11(1 + ^")) - E ^(^) log(x(a)) ] (36) 

(7 a a (7 

The solution can be again found assuming that 2" — 1 quantities a; (a) take the value 
and the state with a°' = 1 for all a takes the value as in formula ( l23l) . 

We can compute as usual and expand to the first order in n, than \E'„ is the maximum 
over 6 of 



nlog(2)| ^2 + ^iog(2)(l - 0) log(l -6)+n log(2)(5 + nlog(2)(l - ^)(1 - log(l - 9)) (37) 

If we take into account fl35|) and the large /3 equalities log(cosh(/3)) = (5 — log(2) and 
log(cosh(/5/i)) = ph — log(2), we can write the energy E as 

E = -^-6h (38) 

while the entropy S is 



S = log(2)^(02 _ ^) ^ i^g(2)(l - ^,) log(l - 5) + log(2)(l - 9,){1 - log(l - 9,)) (39) 
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FIG. 2: Entropy (dashed line) and order parameter Oc (full line) as a function of the dilution 
coefficient 7 at temperature and magnetic field concentration 5 = 0.2. The transition disappears 
and the derivatives of both entropy and 0c are continuous everywhere. 



where 9c is given by the equation 



(l-5)exp(-7^,) = l-^, (40) 

At variance with the magnetic field case, this equation has a single solution 9c. The 
optimum parameter 9c is always positive and the transition disappears since the derivatives 
with respect to 7 of both entropy and 9 are continuous everywhere. In Fig. 2 we plot 
entropy and order parameter 9c as a function of the dilution coefficient 7 at magnetic field 
concentration 6 — 0.2. 

In Fig. 3 we plot the order parameter ^ as a function of the magnetic field concentration 
5 for 5 different values of 7: 7 = 0, 7 = 0.5, 7 = 1, 7 = 1.5, 7 = 2. At S — 0, the order 
parameter 9 vanishes only for 7 < 1 where spontaneous symmetry is broken. 

Finally, in Fig. 4 we plot the entropy as a function of the magnetic field concentration 5 
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FIG. 3: Order parameter 9c as a function of the magnetic field concentration 5. The different 
curves corresponds, starting from below, to dilution coefficient 7 = 0, 7 = 0.5, 7 = 1, 7 = 1.5, 
7 = 2. At 5 = 0, the order parameter 6c vanishes only for 7 < 1. 

for the same 5 different values of 7. The entropy is always smaller then log(2) except when 
both 7 and 5 vanish. Furthermore, the entropy always vanish when 5 = 1 since all spins are 
oriented along the magnetic field. 



7 - DISCUSSION 

In this paper, following Monasson we used magnetization densities as replica order 
parameters. The method is tested against the dilute ferromagnet model even in the case 
in which a randomly located magnetic field is present. We are able to compute exactly the 
ground state energy and entropy. 

When temperature is not vanishing, we only write down the formal solution, while the 
effective one asks for a correct parametrization of the magnetization densities x{a). In this 
case in fact, the problem is not mapped into a simple Potts model. Some work in this 
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direction is in progress. 

The method can be straightforwardly apphed to the dilute spin glass, where we hope to 
find analogous results. Since, in principle, it can be used for all models, without proliferation 

of order parameters, we propose it as a general tool in the replica trick context. 
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FIG. 4: Entropy as a function of the magnetic field concentration S. The different curves corre- 
sponds, starting from above, to dilution coefficient 7 = 0, 7 = 0.5, 7 = 1, 7 = 1.5, 7 = 2. At (5 = 1 
entropy is always vanishing since all spins are oriented along the magnetic field. 
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